Fractional Electromagnetic Waves 
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Abstract: In the present work we consider the electromagnetic wave equa- 
tion in terms of the fractional derivative of the Caputo type. The order of the 
derivative being considered is < 7 < 1. A new parameter a is introduced 
which characterizes the existence of the fractional components in the system. 
We analyze the fractional derivative with respect to time and space, for 7 = 1 
and 7 = 1/2 cases. 

PACS numbers: 03.50.De; 41.20.-q 

The recent interest on the fractional calculus (FC) and in particular in the 
fractional differential equations is stimulated by the applications in various ar- 
eas of physics, chemistry, engineering and bioengineering [l]-[5]. Nevertheless, 
the derivation of such equations from some fundamental laws is not an easy 
matter. The fractional operator reflects intrinsic dissipative processes that are 
sufficiently complicated in nature. Their theoretical relationship with FC is not 
yet fully ascertained. Therefore, it is interesting to analyze a simple physical 
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system and try to understand their complete behavior given by the fractional 
differential equation. 



In this work we will consider the electromagnetic wave equation in terms of 
the fractional derivative of the Caputo type. The solutions to the fractional dif- 
ferential wave equation are given in terms of the Mittag-LefHer function. First, 
we consider the fractional derivative with respect to time and second, the frac- 
tional derivative with respect to space. 

The Maxwell equations for the electromagnetic waves in matter may be 
written as 

f-E = ^-p(r,t), (1) 
V-B = 0, (2) 

VXE = "in' (3) 



c 



^ ^ A-kul^,. . eudE 

VxB = ^ t) + im> (4) 

where p(r, t) and j(r, t) are general, time-dependent distributions of charge den- 
sities and current densities, respectively. In Q D = eE is the dielectric dis- 
placement which is proportional to the electric held E with the electric constant 
permittivity e, and the magnetic held B — /iH, where H is the magnetic field 
intensity and \x the magnetic permeability In the presence of matter the mag- 
netic field intensity H replaces the magnetic induction vector B, in vacuum 
these field quantities are equal to each other H = B. In the case of homoge- 
nous and isotropic medium the parameters e and /i are constants, otherwise are 
vectors. Introducing the potentials, vector A(xi,t) and scalar (f>(xi,t) 

B = V XA, (5) 

E = -- d 4~^ (6) 
c at 

and using the Lorenz gauge condition we obtain the following decoupled differ- 
ential equations for the potentials 

m*o- = (?) 

where ^§ — . v is the velocity of the light in the medium. 

The idea is to write the ordinary differential wave equations (I1I2I3I4P and 
(|7|8[) in the fractional form with respect to t. For this, we propose to change 
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the ordinary time derivative operator by the fractional in the following way 



did? 

> — : — . n - 1 < 7 < n, (9) 

where 7 is an arbitrary parameter which represents the order of the time deriva- 
tive, < 7 < 1, and, (j, is a new parameter representing the fractional time 
components in the system, its dimensionality is the second. In the case 7 = 1 
the expression Q transforms into ordinary time derivative operator 



1 d~< 



= T- (10) 

7=1 dt V ' 



CT 1 -T dt~* 

The following Caputo definition of the fractional derivative will be used [1], 

n-l <7 <ne IV = {1,2,...}, 
where 7 £ iR is the order of the fractional derivative 

f (n) (r) = ^ (12) 

and 

POO 

T(x) = / e-H^dt, (13) 

is the gamma function. 

First, we consider the fractional time derivative. Then, using (|TU1) the 
Maxwell equations ([T]|3]) may be written in terms of the fractional time deriva- 
tives 

V-E = (14) 
V • B = 0, (15) 
VxL -i^^, (16) 

VxB = -Jtj + Jl— __. 17) 



The relations (|5I6I) become 



B = V x 1, (18) 
E = -^-^-V0. (19) 
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Then, applying the Lorentz gauge condition we obtain the corresponding time 
fractional wave equations for the potentials (|7I8|I 



AA -i^d^ A = —r>> (20) 

. . ELL 1 9 2 T , 47T 

In the case, 7 = 1, the equations (I2TJ1) and (|2ip become (J7J and ©. 

If, p = 0, and, j — 0, we have the homogeneous fractional differential equa- 
tions 

= 0. (23) 

We are interested in the analysis of the electromagnetic fields in the medium 
starting from the equations ([22]) and (|23|). We can write the fractional equations 
(j22|) and (123|) in the following compact form 

ax 2 c 2 o^ 1 -^ dt 2 ~< u ' 1 ' 

where z(x,t) represents both A(x, t) and <fi(x,t). We consider a polarized elec- 
tromagnetic wave, then A x (xi,t) = ^4y(xi,t) ^ A z (a;j,t) 7^ 0. The equation 
(|24[) is lineal and a particular solution may be found in the form 

z(x,t) = z e- tkx -u(t), (25) 

where k is the wavevector in the x direction and Zq is a constant. Substituting 
(|2"5J) into we obtain 

+u 2 fc2(T 2(l- 7 ) u(t)=0i (26) 

Redefining 
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w = v 



k 2 a 2 ^ =u 2 a 2 ^, (27) 



where ujq is the fundamental frequency of the electromagnetic wave, the equation 
(fz?6| may be written as 

^ +w yt) = 0. (28) 

The solution of this equation may be found in the form of the power series. The 
solution is 

u(t) = E 2j (-uj 2 t 2 ^, (29) 
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where 



is the Mittag-Lcfflcr function. Substituting the expression (f2U| in (f2l)j) we have 
a particular solution of the equation (f2"4"|) 

js(s, f) = z e- lkx ■E 2 ^(- w 2 t 27 ) . (31) 

In the first case, 7 = 1, the Mittag-Leffler function (|29l) transforms into 
hyperbolic cosines and, from (|27[) u> = lu q . Then 



i?2(— Wq^ 2 ) = cosh/y — cj^ 2 ^ = cosh(ia;oi) = coa(ojQt). (32) 



The expression (|32[) is a periodic function with respect to t. Therefore, in the 
case 7 = 1, the solution to the equation ([24| is 

z(x,t)=Rez e L{uiot - kx \ (33) 

which defines a periodic, with fundamental period To = monochromatic 
wave in the, a;, direction and in time, t. This result is very well known from the 
ordinary electromagnetic waves theory. 

For the second case, 7 = 1/2, the equation ([24)) becomes 
d 2 z(x,t) efi 1 dz(x,t) 



dx 2 c 2 cr <9i 



= 0. (34) 



The solution may be found in the form of (|2"5j) . then we obtain the following 
equation for the function u(t) 

§+^) = 0, (35) 



where, in this case, or = to 2 a, from (|27|) . Solution of the equation (|35|) may be 
obtained in terms of the Mittag-Leffler function (|29|) . In the case, 7 = 1/2, we 
have 

u(t) = Ex{-u 2 t} = er^ 1 . (36) 

The particular solution is 

z(x,t) = z e~ u2t e- lkx , (37) 

For this case the solution is periodic only respect to x and it is not periodic 
with respect to t. The solution represents a plane wave with time decaying 
amplitude. The time in which the amplitude zq decay e times is 

to = —9 = — t~ ■ (38) 
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It is important to note that, 7, is a dimensionless quantity which characterizes 
the order of fractional time derivative while the quantity, a, has dimensions 
of time, and characterizes the presence of fractional time components in the 
medium. However, these two quantities are related as follows 

2 2 

7 = ^o = ^2=t|> 0<a<T o . (39) 

where, T , is the period of the wave, A, is the wavelength and, a = where, 
v, is the velocity of the electromagnetic wave in the, x, direction. 

Taking into account this relation, the solution (|2T)|) may be written as 

u(i)=£2 7 (-7 (1 ~ 7) ? 7 ), (40) 
where, t = -jr, is a dimensionless parameter. 

Now, we will consider the equation (f24|) assuming that the spatial derivative 
is fractional and the time derivative is ordinary. Then, we have the spatial 
fractional equation 

1 d 2S z(x,t) id 2 s(x,t)^ Q 



2(1-5) Q X 2S v 2 dt 2 

where the order of the fractional differential equation is represented by < 6 < 
1, and <j x has length dimension. A particular solution to the equation (|41j) may 
be as follows 

z(x,t) = z e iujt u(x). (42) 
Substituting (gSJ in we obtain 

d ^r + = 0, (43) 

where 

p = ^l a 2(l-8) = k 2 a 2(l-S) (44) 

is the wave- vector in the medium in presence of fractional components, and k 
is the wave vector in the medium without its presence. The wave-vectors are 
equal, k — fc, only in the case, (5 = 1. Solution of the equation f|43[) is given in 
terms of the Mittag-Leffler function 

«w= B -(-* Vi )=ggJri)' <45) 



First case: For the fractional spatial case, when 6=1, from ([44]) and (|45p . we 

have 

E 2 (-k 2 x 2 ) = cosh( v / -fc 2 2; 2 ) = cosh(-ikx) = Re(e" lfca: ). (46) 
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In this case the solution follows from (|4"2"j) 



z(x,t) = Rez e %ut -* kx , (47) 

with, k = k = — , where, k, is the component of the wave- vector in the, x, direc- 
tion and is related with the wavelength by, k — j. The solution (J4TJ) represents 
a periodic, with respect to t and x, monochromatic wave. 

Second case: For the case 6 — 1/2, we have from (|44l) . k 2 = k 2 a x = ^cr, 
and [fc 2 ] = j has dimensions of the inverse of the length. The solution for this 
case, has the form 

u{x) =Ei(- k 2 x) = e-^ x . (48) 



The solution (|421) is written as 

z(x,t) = £ e iut e- px . (49) 

The wave is periodic only with respect to t. The distance at which the amplitude 
zq is reduced e times is 

x ° = T* = ik- (50) 

In this case we have that, 5, is a dimensionless quantity and a x is related to the 
fractional space. These two quantities are related by 



a 



2 



8 = k 2 ai = ^. (51) 
We can use this relation in order to write the equation (1431) as follows 

u{x)=E 2S (-5^Z 2& ), (52) 

where, X — ■ "T j IS 81 dimensionless parameter. It can be seen that the solutions 
(|4"0|) and (|52"|) have the same structure, and in the case, (5 = 7, have the same 
shape. Then, we can plot the function u(s) where s = (x,t), for different values 
of the fractional parameter 7, (see Figure 1). 



Conclusion: In this work we have studied the behavior of the electromag- 
netic waves applying the formalism of the fractional calculus. The order of the 
derivative being considered is < 7 < 1. It showed that for the case where 
7 = 5=1 the solutions represent a periodic, with respect to t and x, monochro- 
matic wave, as it should be. However, if we take 7 = 1/2 the periodicity with 
respect to t is broken and behaves like a wave with time decaying amplitude 
Eq. (|5T|) . On the other hand, when S — 1/2 the periodicity with respect to x is 
broken and behaves like a wave with spatial decaying amplitude Eq. (|4"9l . 
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Figure 1: Graph corresponding to the equations (j40|) and ([52} 
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